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Abstract 

> 

^ The measurement of the 7?^ mass, together with a QCD result for the hyperfine splitting Ehfs = 

2 -^T(i5) ~ ^vb^ allows us to determine the strong coupling constant as at a low energy scale. The 

CN result 

q as{Mr^,s)) = 0.197 ± 0.002|^^e..^ ± 0.002 |^^,_ ± 0.002 |,^^,^ ± 0.006 |,^^ ± 0.005 , 

> 

a,(M^o) = 0.124 ± 0.001 ± 0.001 |^^,_ ± 0.001 |,^^,^ ± 0.003 |,^^ ± 0.002 



ho 

is compatible with the current world average of as reported by the Particle Data Group, and 
shows that the experimental lowest-lying bb hyperfine splitting can be reproduced in terms of a 
perturbative and nonperturbative QCD contribution. 

PACS numbers: 14.40.Pq, 12.38.Qk 
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The observation of the rji, by the BaBar [H |2] and CLEO Collaborations [3] comes after 
three decades of searches of the lightest pseudoscalar bb meson. The spin-singlet state rj^ has 
been detected in the radiative T(3S') — > rj^'j and T(2S') — )■ rj^'y decay modes, studying the 
spectrum of the final photon. The measured mass reported by the BaBar Collaboration is 

Mri, = 9388.9;|^ (stat) ± 2.7 (syst) MeV (1) 

from T(3S') — )■ rjb'y \T], and 

Mri^ = 9394.2+^:^ (stat) ± 2.0 (syst) MeV (2) 

from T{2S) ^ Vbl ^■ 

The signal of rjb in the T(3S') — ?■ 7757 radiative decay has been confirmed by the CLEO 
Collaboration, which quotes [3] 

Mr,, = 9391.8 ± 6.6 (stat) ± 2.0 (syst) MeV . (3) 

The three mass measurements, combining in quadrature the statistic and systematic 
uncertainties, produce the average value [1] 

M^^ = 9390.9 ±2.8 MeV (4) 

which is a remarkable result, since it provides us with a measurement of the hyperfine 
splitting (HPS) of the lowest-lying bb doublet, 

Eh^s = Mrns) - M,, = 69.3 ± 2.8 MeV , (5) 

where we used the Particle Data Group value Mr{is) = 9460.30 ± 0.26 MeV for the mass of 
T(l^) [4J. 

The experimental result (|5| can be compared to the predictions of quark models [5], 
lattice QCD [6], and QCD sum rules [7]. Moreover, it is particularly interesting since it can 
be compared to the result of evaluations based on perturbative QCD with the inclusion of the 
leading nonperturbative contribution, an expression involving fundamental QCD parameters 
such as the strong coupling constant as at a low energy scale. The obtained value of can 
be compared to other determinations, and considered when the average is carried out. In this 
way, one can also investigate if there is room, in the experimental result, for contributions not 
related to QCD, such as that from the mixing effect envisaged in [8] under the assumption 
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of the existence of a light CP-odd pseudoscalar Higgs. The determination of as from the 
result ([5]) is the purpose of the present study. 

As recognized since the early studies of quantum chromodynamics applied to mesons 
comprising heavy quarks a quantitative description in QCD of the QQ bound state 
is possible if the average distance between the quark pair is smaller than the typical QCD 
length scale r ~ 1/ Kqcd'-, the description is given in terms of a perturbative QCD expression 
and nonperturbative corrections. 

The perturbative contribution to HFS comes from diagrams with external heavy quark- 
antiquark lines and the exchange of gluons and light-quark loops. At the leading order in 
the as expansion, such diagrams produce, in the static limit, the Coulombic QQ potential 



with Cp = {N^ — l)/(2A'^c) the eigenvalue of the quadratic Casimir operator of the funda- 
mental representation of the group SU{Nc), Nc being the number of colors. The expression 
of the potential at one- [10] and two-loop [TT] orders has been recently enlarged to three 
loops [12]; the meson masses are obtained as energy levels of a Schrodinger equation. 

At the leading order in a^, the perturbative contribution to the hyperfine splitting is 
proportional to the beauty quark mass and to the fourth power of as{^): 

HFS ^ ' \ J 

In this expression, the dependence of on the renormalization scale /i requires a proper 
choice of this parameter in order to apply the formula to the physical case. A milder n 
dependence can be achieved including higher order corrections. The 0{as) corrections to 
the Ejfpg leading order term have been computed in Refs. [IHl [H], and the result includes a 
logarithmically enhanced as log ag term. Such kinds of terms can be resummed to all orders 
through a renormalization group analysis carried out in the framework of the (potential) 
nonrelativistic QCD effective theory, and indeed in Ref. p5| a Next-to-Leading Log (NLL) 
expression of the hyperfine splitting has been derived which includes a resummation of terms 
of the form a" log"~^ as, with as renormalized in the MS scheme. A discussion can be found 
in [16]. 

In the following we use the formula for E^p^ in [15] (corrected in version 2 of the preprint 



a. 
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in the arXiv), together with the expression of to four loops [T7] : 



1- + 



PI 



(in^L -InL - l) + 



/3o 



+ 



§f_ln3L+-ln2L + 21nL-- 
Pl\ 2 2 



where L = In {^"^ / A^) and /3j are given by [18 
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is the number of active flavors {nj = 4 in the case of the hh system) and ^3 = C(3)- 
The renormalization group improved expression of E^p^ involves the beauty quark mass 
rrtb as an overall factor. It also involves the strong coupling ctg evaluated at a low energy 
renormalization scale /i and at a matching scale m^, as well as on the QCD parameter A. 
The dependence on the two different scales and allows us, through an error analysis, 
to bound the value of A. 

For a heavy QQ pair the nonperturbative contribution to the hyperfine splitting is related 
to the dynamics of the colored quarks in the gluon background. If the size of the quarkonium 
system is smaller than the fluctuations of the background gluon field, this background field 
can be considered homogeneous and constant, and parametrized by gluon condensates. The 
lowest dimensional nonperturbative contribution to the hyperfine splitting of the lightest hh 
doublet has been evaluated in [191 120] and involves the dimension-four gluon condensate. 
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with Pi = (47r)*"'" since the effective Coulombic potential V^j^j = —Cp — has been consid- 
ered with the inclusion of two-loop corrections [13l [21] . The parameters Oi and a2 are given 
by 
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with jE the Euler constant, Ca = Nc the eigenvalue of the quadratic Casimir operator of 
the adjoint representation of SU{Nc), and Tp = |. This contribution must be added to 
the perturbative one, and represents the first term, in the vacuum condensate expansion, of 
a series involving condensates of higher dimension and higher powers of the inverse heavy 
quark mass |22j . 

In the theoretical expression of Ehfs the coupling constant as appears both in the per- 
turbative and nonperturbative terms, and the formula 



E 
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(13) 



involves the factor nib, the renormalization and the matching scales, together with the QCD 
parameter A*^"'-^^^^ the QCD scale in our problem with four active fiavors. We fix the 



Oi 

gluon condensate to the commonly accepted value (0| — G^'^^lO) = (0.012 ± 0.004) GeV^ 



[23] . and we include the uncertainty on the value of nib in the denominator of (10) in the 
uncertainty of the condensate. In the same uncertainty, which is of about 33%, we can 
also include the effect of a possible difference in the scale of in the non perturbative 
contribution with respect to the perturbative one. At odds with other analyses, we keep 
the nonperturbative contribution using its face value, instead of, e.g., fixing it from the 
charmonium hyperfine splitting and then rescaling to the bottomonium case [T5| . 

Proceeding in the numerical analysis, we divide E'^pg and Ehfs by Mx/2, obtaining 
^'hfs ^"^^ Ehfs- In ths QCD expression, this amounts to dividing by the mass rrih defined 
in the IS* bottom quark mass scheme. The change of the scheme induces a higher order 
as correction in the theoretical formula, which is beyond the chosen level of accuracy; in 
Ehfs the remaining rrib dependence is only logarithmic, and encodes the dependence of 
the matching scale. Hence, the procedure is to evolve as from M^o to with four-loop 
accuracy, and below rrib to use as according to the logarithmic precision of the expression for 
Ehfs- To check the uncertainty in this procedure, we have also used full four- loop accuracy 
for as at all scales; the difference in the numerical result is included in the final error budget. 

We define the function 

e(A, /i, m,) = ^^^^^ ^5^^ (14) 

where A = A*^"-''"^^ and search its minima in the space of parameters /x, A and mh- To fix the 

range of A allowed by the comparison between the experimental datum and the theoretical 

expression we proceed in the following way. The minima of ^(A, /i, rri},) determine an implicit 

relation between its variables. Once we have obtained a set of three values (A*,/i*,m^) 

corresponding to the minimum of ^, we fix m?, to the value ml and study ^{A,^,ml) as 

a function of A and yU (the central band depicted in FigjT|. The procedure is repeated 

using in the expression of ^, as the experimental datum, the values E'j^pg — /S.E'^^pg and 

Ehfs + '^■^hfs respectively, obtaining the left and right bands in Fig. [l| hence bounding 

the parameters A and /i to a region of the parameter plane depicted in FigjlJ 

Along the curves of minima, there are ranges of A*^"^^^-* where the dependence on the 

renormalization scale /i is minimized. We bound A^'^^^'^^ in these ranges, imposing, that the 

condition — = is satisfied: the corresponding band is depicted in Fig 1 In this 

Ofl |_ 

way we bound within 1800 MeV < ii < 2200 MeV (slightly larger values than the scale 
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Figure 1: Correlation between the renormalization scale fx and A^"/^^) from the bb hyperfine 
splitting at a fixed matching scale. The central (blue) curve refers to the central value of E^pg, 
the left and right (gray) curves to E'j^p^ — AE'j^p^ and E'^^pg + AE'j^^'pg , respectively, and the 
matching scales are rrih = 4724, 4746 and 4691 MeV, respectively. To minimize the dependence of 

A^-"/ ^' on u, a vertical band is found through the condition = 0. 

Ofi 



/i ~ 1500 MeV chosen in [15] to compute the central value of Ef^pg), and we obtain for 

A{"/=4)^ 



-12 
-13 



MeV 



(15) 



At yU = 2000 MeV the values of A^^'f^^'' where the function ^ vanishes are depicted in Fig|2j In 
this region of the parameter space, the perturbative contribution to the hyperfine splitting 
amounts to E^pg = 65.84 MeV, while the nonperturbative contribution is E^pg = 3.58 
MeV; therefore, for the lowest-lying beauty doublet the splitting is mainly of perturbative 
origin. It is interesting to consider the case where the non perturbative term is forced to be 
zero. In this condition, the experimental value of E'j^pg is reproduced through a larger value 
of the QCD parameter: A("/='^) = 414+}° MeV, so that an uncertainty of 16 MeV can be 
attributed to A*^"^^^^ from the D = A gluon condensate. Moreover, using four-loop accuray 
in as at all scales or the procedure of using as with different accuracies described above 
induces an error of ±20 MeV on A^"-^^^\ quoted as AA*^"'-''"^)| scheme in the final result. 

In the rrih — A parameter plane, a simple correlation is found between the parameters, 
which can be easily understood due to the logarithmic dependence of E^pg on m^/A. The 
uncertainty on A is linked to the variation of the matching scale, which cannot be sensibly 
larger than My/2. A variation of 250 MeV of this scale induces a shift of A^"'-''"^) of about 
40 MeV, an uncertainty dominating the error of A^"-^^^\ The central value of the result for 
^{n/=4) corresponds to rrih = My/ 2. 
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Figure 2: The function ^(A, fi, mi,) at ^ = 2 GeV for the three cases E'^^fs (continuous blue curve) 
and E^Hpg ± AE'^pg (dashed brown curves). The values of nif, are fixed as in Figjl] 

The last source of uncertainty comes from the neglect of (uncalculated) higher order 
contributions to E^fs- The size of these contributions has been estimated considering the 
difference E^p^ — E^^g, with the conclusion that it is about 20.5% of the central value of 
Ehfs [15] • This effect produces an uncertainty, quoted as ho (higher orders), of ±33 MeV to 
7\{»^/=4)_ conservative, we include this uncertainty in the final error, even though some 

higher order effects (for example in the accuracy of as) have been considered separately. 

Following all the steps in the outlined procedure, we obtain a result for A*^"-^^'^^ from the 
experimental bb hyperfine splitting: 



^W=4) ^ 3gg + 12 
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With the value of A^"-^"^) in (16) it is possible to evolve to the T(IS') and to the 
mass scale, implementing the proper matching condition at /i = Mf to include the fifth 
flavor, the beauty [TTl | 



ai'''-'\Mf) = ai'''\Mf) 
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.1 , 11 564731 82043, 2633 , ,x . , 

with k, = - and h = - - ^(-/ - 1). We find 

«.(Mt(i5)) = 0.197±0.002|^^^..^ ±0.002 L,,_±0.002 ±0.006 |,^^ ±0.005 (1^ 

and 



<Mzo) = 0.124±0.001 U^... ±0.001 L,,_±0.001 |,^^,^. 



±0.003 L ±0.002 L . (19) 

om.i. ho ^ ' 
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Figure 3: Measurements of as{Mzo) used in [26] to compute the 2009 world average, togetlier with 
the determination of as{Mzo) obtained in this paper. The continuous vertical Hne corresponds 
to the world average value in [25], and the dashed lines take the error into account: as(M^o) = 
0.1184 lb 0.0007. This result is dominated by the HPQCD determination (indicated as QQ states) 



Equations ( 18 ) and ( 19 ) show the quality of the determination of the strong coupling constant 
from the hh hyperfine splitting: for comparison, the determination of from the ratio 
= r(T — )■ ■ygg)/TiT — )■ ggg) of radiative/hadronic decay widths of T(IS') corresponds 
to a,{Mr) = O.lSAtofil and a,{Mzo) = 0.119+°;°°^ [25]. 



The result in Eq.(19) can be compared to the world average of as- The 2009 average 
computed in [26] is obtained considering, together with the result from the ratio of radia- 
tive/hadronic T decay widths, the determinations of as from r-lepton decays, deep inelastic 
scattering processes (in particular, from nonsinglet structure functions and jet production 
rates), e+e^ processes (event shapes and jet production rates), and from electroweak preci- 
sion fits. Moreover, a determination of the HPQCD Collaboration, based on the analysis of 
the QQ system on the lattice, is included: as(M^o) = 0.1183±0.0008 [27]. As one can see by 
looking at Fig. [sl this last value dominates the present average, as(M^o) = 0.1184 ± 0.0007 



[26] , and compared to this average, the value quoted in (19) is less than 2a higher. The 



result (19) from the bb hyperfine splitting, as{Mzo) = 0.124 ±0.004 (with the error obtained 
by combining in quadrature the various uncertainties), together with the 2009 world average 
of as obtained in [26], slightly increases the value: 

as{Mzo) = 0.1186 ± 0.0007 . (20) 

Our conclusion is that there is the possibility to accommodate the experimental datum 
on the hyperfine splitting of the lowest-lying bb doublet with the QCD, perturbative and 
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nonperturbative, description of it. The resulting value of ccg is compatible with the world 
average within less than 2 standard deviations. The inclusion of the value of ctg determined 
in this paper slightly increases the average. 
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